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Abstract 

We carry out a thorough analysis on a class of cosmological space-times which admit three 
spacelike Kihing vectors of Bianchi class B and contain electromagnetic fields. Using dynamical 
system analysis, we show that a family of electro-vacuum plane-wave solutions of the Einstein- 
Maxwell equations is the stable attractor for expanding universes. Phase dynamics are investigated 
in detail for particular symmetric models. We integrate the system exactly for some special cases to 
Qi confirm the qualitative features. Some of the obtained solutions have not been presented previously 

to the best of our knowledge. Finally, based on those analyses, we discuss the relation between those 
homogeneous models and perturbations of open Friedmann-Lemaitre-Robertson- Walker universes. 
We argue that the electro- vacuum plane-wave modes correspond to a certain long-wavelength limit 
of electromagnetic perturbations. 
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I. INTRODUCTION 



Cosmological magnetic fields have been studied for decades. Various upper bounds liave 



been placed on the strength o 



any primordial galactic and extragalactic magnetic field 
w Faraday rotation measure {l,y], CMB anisotropy 3-5| and primordial nucleosynthesis 
6|, It] although they are inconclusive regarding its existence. Recently, a lower bound of 
B ~ lO^^^G has been claimed by using TeV gamma-ray sources js, 9|. Given the increasing 
quantity and accuracy of the data from astronomical observations, it is important to examine 
various theoretical possibilities concerning the large-scale electromagnetism in curved space- 
time even if the effect is likely to be small. 

It is well known that the energy density of electromagnetic field decays as the inverse 
fourth-power of the scale factor of the universe when it is perturbatively analyzed around 
a flat Friedmann-Lemaitre-Robertson- Walker (FLRW) background due to the conformal 
invariance of the Maxwell equations. Recently, however, it has been pointed out that in open 
FLRW universes, the decay of magnetic energy density can be slower than that of blackb ody 



radiation if you take into account modes with wavelengths above a certain threshold |lO |. 
The source of amplification is the interaction between the scalar curvature of the three-spaces 
and the vector-mode perturbation of electromagnetic field, which would be a second-order 
effect in the perturbation around fiat FLRW. Based on that result, it is interesting to look 
into nonlinear effects of the magneto-curvature coupling on large scales. 

Spatially homogeneous Bianchi cosmologies are a suitable framework to generalize FLRW 
cosmology in the limit of small long- wavelength inhomogeneities. While the model is 
tractable as a system of ordinary differential equations, it is fully nonlinear and compli- 



cated enough to exhibit qualitatively new behaviour. It is known [Uj that only types I, II, 
VIq and VIIq admit pure magnetic fields among all the different models of spatially homo- 
geneous universes and they have been investigated by many authors 12|-|l5|. Although the 
other models have been given less attention in the context of magnetic fields because of the 
lack of observational evidence for large-scale electric fields, we need to investigate them for 
the nonlinear effects described above because only type V and Vllh contain open FLRW 
models as special cases. Indeed, it is this mixing of electric and magnetic mode that gives 
rise to the amplification in the perturbative analysis. We expect, at least intuitively, that 
the isotropic limit of those homogeneous universes should reproduce the behaviour of a long- 



wavelength limit of perturbations around open FLRW. These models belong to a subclass 
of homogeneous cosmological models, the so-called class B of Bianchi cosmologies which is 
of interest from a mathematical point of view as well because of the existence of electro- 
vacuum plane-wave solutions that are generalization of vacuum gravitational plane-waves. 
In the pure-gravitational class B models, the vacuum plane-wave space-times are found to be 
stable attractor solutions for expanding initial conditions [l6|. When a tilted perfect fluid, 
whose peculiar velocity is not aligned to the unit normal of the homogeneous spatial slices, 
is included, it has been shown that they are not necessarily simple attractor solutions and 
can exhibit non- self- similar looplike behaviour at late time [l^]. Since the inclusion of the 
electromagnetic field introduces energy fluxes with respect to the homogenous hypersurfaces 
(as does that of titled fluid), it is natural to ask whether these plane- wave space-times are 
stable in the Einstein-Maxwell system. 

In the present paper, we carry out a thorough analysis of the nonexceptional Bianchi 
class B cosmological models containing a general electromagnetic field and an orthogonal 
(nontilted) perfect fluid with a linear equation of state (we refer to it as a 7-law fluid). 
First, in section II, we reduce the governing equations to a standard form which is suited to 
analyze the stability of various self-similar solutions. In section III, the dynamical system 
analysis suggests that the 2-parameter family of electro-vacuum plane-wave solutions is a 
stable attractor and possesses features very similar to its pure-gravitational counterpart. 
In section IV, we take a closer look on the axisymmetric subcases which still contain open 
FLRW models as a special case. In section V, we integrate the equations analytically in some 
special cases to complement the qualitative analysis and to facilitate comparison with the 
perturbative analysis later. Some of the solutions have not been found before to the best of 
our knowledge. In section VI, we look for a connection between these nonlinear homogeneous 
models and the large-scale limit of linear perturbations and argue that the Bianchi models 
correspond to a certain long-wavelength limit of vector perturbations around an open FLRW 
universe. 



II. MAXWELL EQUATIONS AND THE DYNAMICAL SYSTEM 

The notations and terminology are mostly employed from Wainwright and Ellis 
(chapter 1 for covariant approach to Bianchi cosmologies and chapter 4 for application of 



dynamical systems analysis to cosmology). We adopt the convention where four- vectors are 
presented in bold face and 1-forms come with over-bars. Otherwise all the quantities are 
understood to be real numbers. Latin indices run from to 3 and Greek letters are used to 
label the spatial part (1 to 3) of them. Differentiation by proper (clock-) time t is denoted 
by an over dot. 

n 

Following Ellis and MacCallum [191], we take a group invariant orthonormal frame {ca} 
and their dual 1-forms {w"}. The fundamental variables describing the geometry of the 
space-time are the commutators among the basis vectors 

[ea, efc] = 7"„;,ec 

and their nonvanishing components are prametrised as 

where ea/5-y is the three-dimensional Levi-Civita symbol and Greek indices are raised/lowered 
by 5"^ /dap- A proper time coordinate labehng the homogeneous hypersurfaces is defined by 

d 

and all the variables above can be seen as functions of t since 

ea(7V) = 0. 

As the matter contents of the universe, we take a nontilted 7-law perfect fluid, whose 
energy density is denoted by p, and a source-free electromagnetic field. It is natural, though 
not always necessary, to demand that all the components of the electromagnetic field in this 



frame depend only on t. In defining the 1 -|- 3 s 
convention of Misner, Thorne and Wheeler 



Dlit of field strength tensor Fab, we follow the 



20(; that is 



The source-free Maxwell equations 



dd = d*d = 
4 



can be written in terms of the components in the orthonormal frame with a help of the 
relation 

The result reads as follows: 

= -2HHa + <y^pH^ + e^p^H^W + n^pE^ + e^^p^ E\ 
K = -2HE^ + aapE^ + e^p^E^n^ - n^pH^ - e^p^ H\ (1) 

The isometry group G3 of Bianchi class B admits an Abelian subgroup G2 and we choose 
62 and 63 to be tangent to the orbits of the 6*2- Then the constraint equations in ([1]) imply 
El = Hi = 0. In the Binachi class B space-time containing this source-free electromagnetic 
field and the nontilted perfect fluid, it turns out that the G2 acts orthogonally transitively 

led nonexceptional case. 

Thus we can use the equations given in the section 1.6.3 of Ref. 
The source terms are taken to be 



18l setting di = iii = 0. 



Id = p + 37r+, 

P = (7 - 1)P + 
1 

^\i.^2 + Hi - El - Hi) -{E2E3 + H2H2, 



T,^= [El + El + Hl + Hl), 



T^AB = 

qi = E2H3 — E3H2 



iE2Es + H2Hs) \{El + HI - El - Hi) 



We assume 7 is constant and satisfies < 7 < 2. 

Now we follow a standard procedure established by Wainwright and his collaborators 



2l| to reduce the equations into a form suited for qualitative analysis. In doing so, we 
carry out the 1 + 1 + 2 split of the space-time so that the gauge freedom about the 1-axis 
is manifest. We classify the expansion normalised variables according to their behaviour 
under a rotation around the 1-axis. 



Scalar s (spin-0) : 



3if2' H' " 3if2- 
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Tensors (spin- 2) : 





0"_ 




n_ 


n_ = 




CTx 


iVx 


nx 





7^33 



Inhomogeneous (rotation angle itself): 

The variables 11+ _x and S are quadratures of electromagnetic field components and their 
evolution equations can be obtained from the Maxwell's equations ([1]). We introduce a new 
time coordinate r and deceleration parameter q by 

and use primes to denote derivatives with respect to f. We derive the following equations: 



Einstein equations: 



q = 2(1:1 + + s'x) + ^(37 - 2)n + n+, (2) 
1 = n + + + + + Nl + Nl + n+, (3) 

= H + 2(S+A + S_Arx -SxA^-), (4) 
S', = (g-2)S+-2(iV!+iV2)-n+, 

S'_ = (g-2)S_ + 2Sxi?-2(A^+A^_ - AA^x) -n_, (5) 
= (g-2)Sx -2i:_R-2{N+Nx + AiV_) -Hx- 

Jacobi identities: 

= (g + 2S+)Ar+ + 6(S_iV_ + SxArj, 
A^: = (g + 2S+)iV_ +2(S_iV+ + Ar^i?), 
iV^ = (g + 2S+)Arx +2(SxAr+-iV_i?), 
A' = (g + 2S+)A 

C/ass i? first integral: 

A" = h{Nl-3{Nl + Nl)}. (7) 
Energy conservation for the fluid: 

n' = {2q - 37 + 2)n. (8) 



(6) 



(9) 



Maxwell's equations: 

n; = 2(g - 1 + s+)n+ + 2(s_n_ + SxHx + ae), 

E' = 2(g - 1 + + 2{AU+ + AT^n. - N^U^), 

n'_ = 2(g - 1 + s+)n_ + 6S_n+ - eiv^s + 2i?nx, 
n'^ = 2(g - 1 + s+)nx + 6Sxn+ + 6N^e - 2i?n_. 

There is an algebraic constraint among the quadratures: 

ul = E' + ^{n'_ + nl) (10) 

Here we comment on the degrees of freedom. The dynamical system is not closed because 
there is a remaining gauge freedom as we mentioned above. In order to decide the frame to 
fix the gauge, we face a similar problem as was seen in the analysis of tilted Bianchi class 
B models by Coley and Hervik 171]. Namely, that there is not a simple choice of frame 
for which the governing equations become regular everywhere. Which choice is convenient 
depends on the equilibrium points or invariant sets under consideration and we will make use 
of several different choices in the present article. Note that construction of scalar variables 
such as S_A^_ + SxA^x would lead to too many constraints as was the case in the Ref. [17. 



Whichever choice we make, however, the dimension of the dynamical system is seven. The 
counting goes as follows; we take S+ _^x, A^+,-,x and n+^„_x as the fundamental variables; 
A, E and Q are determined by (171), ( fTOl) and the Hamiltonian constraint ([3]) respectively; R 
will be determined by fixing the gauge, which also imposes an additional constraint among 
the spin-2 variables; finally taking into account the momentum constraint (jl]), we are left 
with nine variables and two constraints. 

III. STABILITY OF PLANE- WAVE SOLUTIONS 

The purpose of this section is to extend the result obtained by Hewitt and Wainwright 

n 

[Iq about the stability of vacuum gravitational waves to Einstein-Maxwell equations. We 
follow the usual procedure of dynamical system analysis. 

A. Invariant sets 



To facilitate the understanding of the dynamical system, it is worth sorting out possible 
invariant sets, which are defined to be subsets consisting of trajectories specified by certain 
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restrictions which form dynamical systems by themselves. Note that we may assume A> 
without loss of generality because of the reflection symmetry about the 2-3 plane. We use 
the integration constant /i = 1//?, of ([7]) to classify the class B models. 

TABLE I. Electromagnetic Bianchi invariant sets. The entries in the upper half are the general 



class B models for the 


source-free electromagnetic field. 




Notation 


Restrictions 


Dimension 




h<0,A>0 


7 


M(Vllh) 


h> 0,A> 0,N+ > 


7 


M(IV) 


h = 0,A> 0,N+ > 0,iV2 = 3{N^+N^) 


7 


M(V) 


h = 0,A>0,N+ = N^=Ny,=0 


5 


M(II) 


A = 0,Nl = 3{N^+Nl) 


6 


M(I) 


A = N+ = N- = Nyc = ^ = 


4 



Table [T] lists the electromagnetic Bianchi invariant sets described by the equations given 
in the previous section. The restriction on the sign of in type Vllh and IV comes from 
the fact that = by itself forces the system to fall into an invariant set and there is 
a discrete symmetry interchanging 2- and 3-directions. We refrain from further discussion 
for disconnected components as it is subtle before the gauge is fixed. M(V) forms a part of 
the boundary of M(1V). M(ll) is the class A boundary for all Vlh, Vll/i and IV while it 



is not the most general e 
investigated by LeBlanc 



ectromagnetic type II. Its four dimensional subset and M(T) were 



14 



TABLE II. Electromagnetic Bianchi class B invariant sets with higher symmetry. 

Notation Restrictions Dimension 

SM{Ylh) h<0,A>0,N+ = 4 

SM{Yllh) /i > 0,S_,x = A^-,x = n_,x = 0,^ > 0,Af+ > 2 

5M(V) /i = o,S_,x =A^±,x =n_,x =0,^ > 2 



Some symmetric invariant sets are shown in Table [TTl 'S'M(VIft,) has four dimensions 
and can support non-null Maxwell fields. SM{Vllh) and SM{V) have an identical two- 
dimensional structure as dynamical systems and the space-time is locally rotationally sym- 
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metric (LRS) 22|. They will play a central role in the comparison to perturbations around 
open FLRW models because we can fully derive the global behaviour of the orbits and 
integrate them explicitly for some interesting cases. 

Aside from the invariant sets listed above, there are obvious electro-vacuum subsets of 
them obtained hj Q = 0. Here we are not concerned with pure-gravitational orthogonal 
class B models (i.e. 11+ = 0) as they were studied by Hewitt and Wainwright [l^. 



B. Equilibrium points 

In dynamical systems analysis, equilibrium points, which are defined as time-independent 
solutions of the system, play an important role. Here we focus on two of them which are 
important concerning the past- and future-asymptotic behaviour of the system. The full list 
of the equilibrium points in the present setup is given in table IIIII 

Kasner equilibrium points /C: Those are a one-parameter family of equilibrium points 
corresponding to the well-known Kasner vacuum solutions. They have always been the 
key to understand past-asymptotic behaviour in Bianchi cosmologies and we expect 
the same is true here. To carry out the stability analysis, it is convenient to fix the 
gauge by setting Sx =0. i? is determined by 

E_i? = -A^+A^x - ^A^- - -Hx. 
In the resulting system of equations, /C is characterised as follows: 

S+ = costp S_ = sin?/' — 7T < < IT 
A = A^±,x =U+ = Q = 0. 

This can be pictured as a unit circle on plane with ip measuring the angle. 

Since all the internal curvatures (A, A^j. x) and electromagnetic quadratures are zero, 
the constraints ([7j) and flTUl) become degenerate. Therefore we have nine eigenvalues 
as follows: 



Ai = A2 = -3(7-2) A3,4 = 2(l + cos^±y3sinV') 
-^5,6,7 = 2(1 + cosip) A8,9 = 1 + cosijj ± y/Ssmip. 
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TABLE III. Equilibrium points in the electromagnetic Bianchi class B. Some of the symbols rep- 
resent several distinct points in the state space because they lie in disconnected boundaries with a 
same structure. For the notation of the invariant sets, see Hewitt and Wainwright {l6^ 



Symbol 


Invariant set 


Self-similar solution 


Restrictions 


P(I) 


B{1) 


Flat FLRW 




P(II) 


B{ll) 


Collins-Stewart 


i<7<2 




B{Ylh) 


Collins 


^>3^'i<7<2 


PMi(I) 


M(I) 


Jacobs magnetic I 


|<7<2 


PM2(I) 


M(I) 


LeBlanc 


f <7<2 


PAf(II) 


M(II) 


Dunn-Tupper 


|<7<2 


/C 


B{1) 


Kasner vacuum 




J 


B{1) 

M{\lh) 


Jacobs stiff fluid 


7 = 2 
/i< 


VM{h) 


M{\llh) 
M(IV) 


Electro-vacuum plane-waves 


/i > 
/i = 


M{h) 


SM(VII,,) 
SM{\) 


Milne universe 


/i > 
/i = 




SM{\llh) 

SMiy) 


Open FLRW 


/i > „ 

'7 = 1 

/i = 



It is difficult to decide which of them represent the true dynamical degrees of freedom. 
However, as we can see, the arc specified by — | < ip < which we shall call /C, has 
positive eigenvalues for all the directions aside from the zero around the circumference. 
The rest of K, except for a finite number of points have a stable manifold of at least one 
dimension represented by either Rosen orbits, Taub orbits, or rotating Kasner orbits 
studied by LeBlanc for magnetic type I [l^ and type II 15|. This past-stability of 



the arc K, also suggests that the Mixmaster oscillation will not occur in this class of 
models. We denote the points ip = Qhy Qi and ip = ti hy Ti. They will appear in the 
LRS invariant sets considered later. 

Electromagnetic plane- wave equilibrium points VM^: They are essentially the only 

10 



nontrivial electromagnetic class B equilibrium points and of interest since the pure- 
gravitational counterpart (a subset of the electromagnetic ones) is a stable attractor 



of the expanding class B space-times 16|. The electromagnetic field is null and the 



corresponding self-similar solutions were discussed by Harvey et.al. 



Skea 



23|, Araujo and 



24| , and Hervik |25j in different contexts. Here we use the gauge freedom to take 



iVx = 0, which implies 

They are expressed as a plane of equilibria, namely a two-parameter family of points 
characterised by constants r and s. 

E+ = -r Ex = -N_ = s E_ = Nl = h{l - + 
A = l-r H+ = H = 2r(l -r) - H_ = Hx = = 

0<r<l s^<r(l-r). 

For VI/i, we have an additional restriction 

Thus the area of the plane becomes smaller for greater h. When = r(l — r), 
they correspond to the pure-gravitational plane-wave space-times. The eigenvalues 
are given by 

Ai = 4r - 37 + 2 A2 = A3 = 



A4,5,6,7 = 2(r - 1) ± y -2(5iV2 _^ ^^1+) ± QN+^ Nl + 6H+ 

where the square root always gives pure imaginary values. The two zero eigenvalues 
reflect the fact that the equilibrium points form a two-dimensional subset. For type 
VII/i and IV with | < 7 < 2, they always have a stable subset specified by r < 
1(87 — 2). For type VI/^, this condition cannot be satisfied ii h > 3?^, in which case 
P(Vl/i) becomes the attractor of the dynamical system. These features are very much 
analogous to those found for the pure-gravitational case. Note that the direction of 
Ai corresponds to the perturbation in Q. Therefore, VM{h) is a local sink in the 
electro-vacuum {Q = 0) class B models. 



11 



C. Monotone function 



The condition 17 = gives rise to six- dimensional boundaries of the full electromagnetic 
class B models. Analogous to the vacuum Einstein models we can show that Z = 
(1 + — is monotonic: 

Z' = 2{q - 2)Z - 3(7 - 2)(1 + 

which implies \Z\ is monotone decreasing for 17 = since g < 2 in the interior of M(IV), 
M(VI/i) and M(VIIft). In fact, Z itself is monotone decreasing since Z > 0. To see this, 
first note that (ED ± g]) yield 

1 - = (S+ ± + (S_ ± N^)^ + (Ex T N_)^ + n+ ± S. 

Taking into account (fTOj) . |H| < n_|_ and therefore |S+ ± A| < 1. Then it follows that 

Z = (1 + S+ + A){1 + S+ - A) > 0. 

Thus, we obtain the following result, which is a direct generalisation of the proposition 5.1 



in Ref . 



16 



Proposition III.l. 

For the vacuum Einstein- Maxwell Bianchi class B models, i.e. M{TV)\^_^, M(Vlh)\^_^ and 
M(VII/j)|^_p except for a set of measure zero, the past attractor is IC and the future attractor 
is VM{h). 

Proof. To obtain an everywhere well-defined bounded dynamical system, we partially fix the 
gauge by setting i? = 0. Then we can apply LaSalle's invariance principle and see that the 
past- and future-asymptotic sets are contained in the subset Z' = 0. Since Z is monotone 
decreasing, the past and future are specified by g = 2 and Z = respectively. It is easy to 
see that q = 2 uniquely characterises the Kasner circle IC and the condition (1 + S_|_)^ = is 
consistent with the plane-wave solutions VM{h) only. Note that these characterisations are 
gauge-invariant. Therefore, we conclude that the future attractor is VM{h) as it is always 
stable in the vacuum invariant sets. The past attractor is contained in /C. We can specify it 
to be IC since for all the points in /C\/C except Ti, the only zero-eigenvalue is the one in the 
direction along the circumference. Then the orbits emanating from /C\/C lie in its unstable 
manifold, therefore are at most of measure zero (e.g. see Theorem 4.1 in Aulbach 26|]). Ti 
cannot be a past attractor as Ti G VM[h). □ 
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IV. LRS MODELS 



The LRS assumption leads to a two-dimensional dynamical system SM{Y) or SM{Vllh), 
both of which contain the open FLRW model. Only the null Maxwell field is consistent with 
the geometry. We visualise the dynamics to facilitate the comparison to perturbations 
around FLRW later. The equations are given by 

EV = {q- 2)E+ - n+, 
A' = {q + 2E+)v4, 
q = i(37 - 2)(1 - A') - i(37 - 6)E^ - ^(87 - 4)n+, 

= 4a2e^, 
^] = 1 - ^2 - E^ - n+, 
= VIa. 

The fourth equation was obtained by combining the momentum constraint with (11 01) . The 
physical region in the yl-E+ plane is defined by 

o< A<i, n+ > 0, Q>0 

and its boundaries. There are two disconnected invariant sets E+ > and E+ < separated 
by the open FLRW orbit E+ = 114, = 0. Within each of their interior, we can solve the 
quadratic constraint as 11+ = ±2AE+. Substituting this into the Hamiltonian constraint 
and using > yield 

^ + E+<1, A-E+<1. 

The projection of VM (h) onto the LRS invariant sets is the line A — E_|_ = L At the ends of 
the line are located M{h) and Ti. They have one zero-eigenvalue along the line and another 
4r — 87 + 2 representing the perfect fluid mode. Thus the points with r > |(37 — 2) are 
local sources and those with r < ^{3^ — 2) are local sinks; increasing 7 means increasing the 
unstable part in the line A — E_|_ = 1. The other source and sinks are summarised in table 

m 

We can flnd monotonic functions for both of those invariant sets E+ ^ when 7 > §• 

A' = {q + 2E+)A > for E+ > 0, 
(aA - /3E+)' = [a^^ - /32E+ + (aA + /3E+)2] [1 - (A - E+)] > for E+ < 0, 
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TABLE IV. The list of local sources and sinks in the LRS models. The system shows different 
behaviours according to the direction of Poynting vector. 



Invariant set 


S+ > 


S+ < 


7 


[0,1) 


2 
3 


(1,2) 


2 


[0,1) 


2 
3 


(1,2) 


2 


Local sources 


Qi 


J 


VM{h) 


J 


Local sinks 


P{1) 




Mill) 


P{1) 




VM{h) 



a ^-(37 -2), /3^-(37-6). 

Thus, by the LaSalle's invariance principle, the local sources are the past attractors and the 
local sinks are the future attractors. The phase portrait of a representative case (7 = 1) is 
given in figure [TJ 

The case S+ < with | < 7 < 2 is worth elucidating. Since two distinct orbits never 
meet each other, we infer that an orbit started from a more anisotropic state (larger |S+|) 
than the other ends up in a less anisotropic future-asymptotic state. For those models, both 
past and future attractors have non-zero 11+ but they isotropise and become close to flat at 
intermediate times. 




FIG. 1. Phase portrait of type V and VII^j LRS models for dust (7 = 1). The diagonal dotted line 
represents the projection of equilibrium points VM^^ onto the LRS invariant set. 
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V. EXACT SOLUTIONS FOR LRS SPACE-TIMES 



We are going to present some exact solutions, including new ones, for the dynamical 
system considered in the previous section, which should reproduce the asymptotic behavior 
described there. 

In the parametrisation introduced by Siklos the line element for LRS type V and 
Vllh is given by 

ds^ = -a-Hr^ + a-^dx^ + e^^^"^') {dy"" + dz^) a > 0. 
Defining new variables 

a = - — , = :r' (11) 

a dr dr 

the Einstein equations reduce to 

^ _ 2^ - 2a/3 - 2/32 = 1(37 - 2)0"^ + 3a-\+, (12) 
dr dr 2 

2{a + 13) = a~^qi, (13) 

-2al3 -2 = 1(2- -f)a-^p + 3a-^TT+, (14) 
dr 2 

f + 2/3^ -2 = 1(2- 7)a-V (15) 
Combining eqns. ( !T2|) + (|T^ + 2 x f lTSj) gives the generalized Friedmann equation 

(3^ -2al3 -3 = a-^{p + 3n+). (16) 

The energy density of the electromagnetic field pem is given by pem = 37r+. Note that 
Pem — I1 because only null-field is consistent with the LRS symmetry. We term qi = p^m be 
left-handed and qi = —pem be right-handed. 

A. Left-handed solutions 
Eelctro- vacuum: 

X^ds' = e^'^i-dr' + dx') + e^^^~^\dy' + dz^) 

—00 < r < cxD p > 1 
Alternative metric : ds^ = -df + {vtfdx^ + {vt)^ e'"^"" [dy"^ + dz^) 
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This is just LRS specialisation of VM{h) with s — 0. The parameter u is related to r 
in the section 3 by 



u-1 
iy + 2' 



(17) 



Perfect fluid 7 = |: 



2rT 

W = , ^ i-dr^ + dx") + e2(— )(1 + e2(i-'-)-)((iy2 + dz") 

Vl + e2(i-0r ^ ^ ^ A y 



— oo < r < oo r > 1 



Asymptotic behaviour : 



T^+oo:ds^^ -dt" + (ri)2 [dx" + e-'^(ciy2 + ^^2^] 



Perfect fluid 7 = 2: 



A'ds' = e""-^ sinh^+'" 2T(-dT' + dx") + sinh 2Te-''%dy' + d^') 



0<r<oo r>0 



Asymptotic behaviour : 



T ^0:ds^ 



-de + t'-^dx^ + t^e-2'^"(%2 + (i;^') 



T ^ 00 : ds^ = -dt^ + [(4r + l)t]^ da;^ + [(4r + l)t]^ e-''''{dy' + (i;^') 



Perfect fluid 7 = |: 



-,4rr] 



r — tanhr/ 



(r^ — l)'^drf dx"^ 
4(r cosh Tj — sinh 77)2 cosh^ r] 



r — tanh rj 



—00 < 77 < 00 r > 1 
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The time coordinate rj is related to r by 



— f — v'r^ — 1 tanh ^2\/r^ — It;^ . 



Asymptotic behaviour : 



rj ^ +00 -.ds^^ -dt^ + (^^^^^--p-t^ ' dx^ + (^^^^^--^t^ e-2-(ciy2 + dz^) 



B. Right-handed solutions 
Electro- vacuum: 

X^ds^ = e^^ sinh-5 2r(-dr^ + dx^) + sinh 2Te-^^(dy^ + dz^) 

< T < 00 

Asymptotic behaviour : 



T 00 -.dx"^ = -dt^ + [rfx^ + e'^''{dy^ + d^^)] 



Perfect fluid 7 - ^ 



3" 

A^cis^ = ^ (-dr^ + dx^) + e2(— )(1 - e-2('-+i)-)(dy2 + dz^^ 

^1 _ e-2{r-+l)r ^ ^ ^ 

0<T<oo r>l 

Asymptotic behaviour : 

T^Q:ds^ = -dt" + r idx' + iie-2^^((iy2 + dz^) 

T ^ 00 : ds^ = -di^ + ^rtf [dx'' + e-2^(dy2 + dz'')] 
17 



Perfect fluid 7 = 2: 



X'ds^ = e^'^ sinh^-"'- 2T{-dr^ + dx^) + sinh 2re-''%dy' + dz') 



0<r<oo 0<r< 



Asymptotic behaviour : 



T ^0:ds^ 



-de + t'-^dx"" + i5^e-2-(dy2 + 



T oo : ds^ = -di^ + t^ [dx"^ + e'^'^idy'^ + dz"^)] 



Perfect fluid 7 = |: 



X^ds^ = eS'-''A(7?)-2S(7?) {-B{r]fdrf + dx^) + e^''^-^^A{ri)B{ri) {dy^ + dz^) 



. , , sin (2\/l — r'^ri) , ^ 



1 



VT^-^cot (2-\/l — r^v) ~ f 



0<V<Voo B{r]oo) ^ = r > 
Here we used the convention tliat 



sm ex 



= < 



sin x e — 1 

X e = 
sinh X e = i 



and the corresponding expressions for cos x. The original time coordinate r is given 

by 



^ = — r + \/l — cot ^2Vl — r^T^j 



Asymptotic behaviour 



77 ^ : ds' = -dr + t-^dx' + tse-"^{dy^ + dz') 



r] ^ r]oo ■■ ds^ ^ -dt^ + t^ [dx^ + e'^'^idy^ + dz^)] 
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The vacuum and 7 = 2 solutions were discovered by Ftaclas and Cohen 28j; Roy and 



Singh 291] gave the right-handed radiation solution. The other solutions appear to be new. 
It can be observed that all the left-handed solutions except 7 = | are asymptotically plane- 
waves in the future and the decay rate of electromagnetic energy density is given by t-^. In 
terms of the average scale factor I, it can be expressed as pem oc I 2+1- using the parameter 
in (fT7|) . When they are close to Milne, namely, in the isotropic limit u ^ 1, the decay is 
as slow as /~^. On the other hand, right-handed solutions which are attracted towards the 
Milne universe have a faster decaying rate of The difference between the helicity is 
clearer when 7 = §• Looking at the left- and right-orbit that have a same attractor in J^{h), 
their decaying rate always differs by ^ with r > 1. Generally speaking, as far as future 
asymptotic behaviour is concerned, left-handed solutions have a slower decaying rate than 
right-handed ones. 



VI. COMPARISON WITH THE PERTURBATIVE ANALYSIS 

Having shown the general stability of plane-waves and derived some detailed behaviour 
of the orbits in the LRS models, we now compare them to the results obtained in the pertur- 
bation around open FLRW. The type V and Vllh models are anisotropic generalisations of 
open FLRW and their isotropic limit should exhibit some features seen in the perturbation 
of long wavelengths. We will see that the LRS modes appear to be generic and correspond 
to a super-adiabitc mode discussed by Barrow and Tsagas flo| . 

First, let us review the vector modes of electromagnetic perturbation around an open 
FLRW background. Following the convention of Barrow et.al. jsol, the linearised Maxwell's 
equations in an arbitrary background frame are given by 

Ba = -2HBa - cm\Ea, 
Ea = -2HE, + curlE,. 

For the comparison with the analysis in the present article, we take the frame to be or- 
thonormal and use the time coordinate r introduced in the previous section, which is 
now identified with the conformal time of the Friedmann background. Note that the curl 
terms include the contribution from the isotropic background spatial curvature. Assuming 
Ba{T,x^) = Bk{T)Y^{xi) where a vectorial eigenfunction of the Laplace-Beltrami operator 
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satisfies 



the solutions are given by 



(A + ^)Y^ = dWY' = Y^ = 0, 



B,ir) = ^ev^^ + ^e-^^ (18) 



where a and /3 are integration constants. The magnetic field follows the standard adiabatic 
decay law Ba oc for short-wavelength modes k"^ > 2, while the decay rate is reduced 



on large scales where the wave-numbers are smaller than the threshold /c^ = 2 10|]. For 
the vector perturbation, the magnetic field is always accompanied by electric field, which 
must be generated through the curl term. Substituting the magnetic solution f|T8|) into the 
Maxwell's equations and using the vector identity on hyperbolic space 

curl curl = — (a + ^ j + grad div. 



we obtain 



E,iT,x,) = ( ^e^^ - L'^- ] curl Y„K 



The apparent /-dependence of the electric field is misleading since the curl operator here 
includes another factor of The physical amplitude of the electric field in the orthonormal 
frame behaves it should because of the symmetry between source-free 

electric and magnetic field. In terms of the energy density pem = \{,BaB°- + EaE"-), the 
result is summarised as follows: There are two distinct behaviours above and below the 
threshold fc^ = 2. For k'^ > 2, we just have the adiabatic decaying law pem oc But 



for longer wavelengths, the decay rate is modified to pem oc /-4g±2V2^^r gg^Q]-^ Qf which 
corresponds to growing (/3 = 0) and decaying [a = 0) mode. If the background is Milne, 
then / oc and we have 



p_ oc /-4±2vCT. (19) 

It later turns out to be helpful to compute the Poynting vectors associated with the 
solutions derived above. They are given by 

This means that for each given Fj^, the growing and decaying modes have exactly antiparallel 
Poynting vectors, which is reminiscent of the LRS models where the system shows different 
asymptotic behaviours according to the handedness of the field. 
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Let us now consider the electromagnetic Bianchi models near isotropy. For the limit of 
isotropy, we assume |S-|-^x| <^ 1 in the dynamical system (|2]) - ffTOj) . This implies A^_,x,n-|-^x 
and S are all first order as well. The linearized electromagnetic sector obeys 

2(g - l)n+ + 2AS (20) 

2(g-l)S + 2An+ (21) 

2(g-l)n_ (22) 

2(g-l)nx. (23) 

First of all, we note that the evolution equations for 11+ and S are closed among themselves 
and do not receive any contribution from n_ or IIx , which is exactly the feature of null elec- 
tromagnetic field appearing in the LRS models. That is, any electromagnetic perturbation 
that must have nonzero energy density 11+ inevitably excites the energy flux perturbation S 
as well, while n_ and IIx are independent of 11+ as the quadratic constraint ( ITOj) is degener- 
ate at linear order. The latter two modes are also subject to the remaining gauge freedom. 
Therefore their physical significance is questionable. Fortunately, we can consider the LRS 
modes n+-S separately as far as the isotropic limit is concerned. Moreover, the evolution of 
n_ X is trivial anyway since the equations fl22|) and fl23|) solve regardless of the background 
as 

n_,x oc H-^-^ 

where H is the background Hubble expansion rate. This is just the adiabatic decaying law 
of radiation. From now on, therefore, we focus our attention on the LRS models for which 
we have some exact solutions and full dynamical picture. 

For I < 7 < 2, the generic behaviour of the LRS models is to start from an anisotropic 
singularity dominated by extrinsic curvature, possibly becoming almost isotropic for an 
intermediate time interval and attracted towards electro-vacuum, which is driven by the 
intrinsic curvature synchronised with the null Maxwell field in the future. The handedness 
of the electromagnetic wavelet, namely the direction of the Poynting vector, is important in 
deciding the asymptotic states. In left-handed cases, both past and future dynamics receive 
significant contributions from the Maxwell field appearing as plane-wave space-times while 
right-handed fields only affect the intermediate evolution. 

In the cosmological context, we are interested in the future asymptotic behaviour of 
the electromagnetic energy density pem in the isotropic limit. That is, when an orbit is 
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n'_ = 
n'x = 




FIG. 2. Phase portrait of type V and VII^, LRS models for radiation (7 = |). All the orbits 
achieve their most isotropic state at late time. 

attracted towards Milne point M{h). This is also the situation where the magnitude of 
the shear variable takes its minimum in most cases: i.e. all the right-handed models 
and left-handed models with 7 < | (FIG. [2]). We can solve the perturbed LRS Maxwell's 
equations fl20|) and fl2Tl) with q = 0, A = 1 and obtain the two independent solutions 

n+ oc oc + const.. 

They translate into the decaying mode pem and the growing mode pem Looking 

at the exact solutions for 7 = | and 7 = 2, we can easily see that the decaying (growing) 
mode represents the right-handed (left-handed) orbits. In fact, this is the generic feature of 
asymptotic behaviour of the LRS models in isotropic limit. Near isotropy, the decaying rate 
of electromagnetic energy density is decided by the handedness of the field. Recalling the 
perturbative result for open FLRW (fTOl) . the long- wavelength inhomogeneous perturbation 
of fc^ = 1 shares the decay rates and the relation to the handedness with the LRS modes 
studied here. This appears to imply that the super-adiabitc mode of /c^ = 1 in open FLRW 
generalizes to the left-handed LRS orbit attracted towards the close-to-isotropy plane-wave 
solutions. 

We can see another isotropic asymptotic behaviour for 7 = | (FIG. [3]). The future- 
asymptotic state of the Bianchi model is J-'{h) and therefore we should compare it with the 
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FIG. 3. Phase portrait of type V and VII/i LRS models for 7 = |. The orbits are attracted towards 
the hne of equilibrium points J-{h) indicated by the vertical dotted lines. 

perturbation around the one-parameter family of open FLRW backgrounds J^{h) instead of 
the Milne universe. Looking at the LRS exact solutions, we see that the electromagnetic 
energy density goes as for right-handed solutions and for left-handed ones with 

r > 1. Here r parametrizes different members within the family of J^{h) as well as different 
obits. On the other hand, the inhomogeneous perturbation around the FLRW backgrounds 
J^{h), labeled by the same parameter r, gives the growing mode 

Pem OC OC / r 



and the corresponding decaying mode is 



Pem OC OC / - . 

Again, we can identify the k"^ = 1 growing mode with left-handed asymptotic behaviour and 
the decaying mode with the right-handed one. The presence of the isotropic perfect fluid 
smears out the effect of magneto-curvature coupling. 

It should be mentioned that these arguments are far from conclusive since the orbits 
that are future-asymptotic to Milne are of measure zero in the general electromagnetic class 
B models. On the other hand, there is no LRS left-handed orbit which approaches Milne 

23 



and we can only have approximately isotropic final states there while all the right-handed 
orbits approach M{h). In this sense, the above comparison is asymmetric between left 
and right, in turn growing and decaying modes. However, it simply indicates that the 
nonlinear interaction can lead to a quite different behaviour from linear perturbations and 
the dominant effect their seems to be the electromagnetic plane-wave mode. It is reasonable 
that this nonlinear plane-wave corresponds to the growing mode of the perturbation while 
the decaying mode can only be seen in the exceptional orbits which are asymptotic to Milne. 



VII. CONCLUDING REMARKS 

We have investigated a class of spatially homogeneous Einstein-Maxwell space-times and 
described possible asymptotic behaviours by a dynamical systems analysis. The results are 
analogous to pure-gravitataional models, with the electromagnetic field acting as a kind of 
bridge between extrinsic and intrinsic-curvature dominated regimes. The extended electro- 
magnetic plane-waves are stable attractors of the system. In the LRS models, we derived 
more detailed features of the dynamics by looking at some exact solutions. The handedness 
of the null field plays a crucial role in the magneto-curvature coupling in the restricted class 
of models. These LRS null Maxwell modes appear to generalise the electromagnetic vector 
perturbations around open FLRW with the wave- vector k = 1. 

The Bianchi models provide another example of the close relation between gravity and 
electromagnetism. The dynamics here are surprisingly simple considering the dimensionality 
of the system. It is interesting to note that Maxwell fields can dominate over perfect fluids, 
for example dust at late times in a long- wavelength limit, even though the simple adiabatic 
decaying law of electromagnetic energy density suggests otherwise. 

As to the connection to the perturbations, this analysis shows that we might not nec- 
essarily be able to ignore the vector mode with the wave-number smaller than a/2 in open 
FLRW models as used to be done (e.g. Ref. (sil). It was already argued for scalar modes 
that we should take into account the supercurvature mode k < 1 when we are concerned 
with a random distribution of perturbations over the sky, even though any causal perturba- 
tion could be described by the subcurvature modes which span the basis of square-integrable 
functions over hyperbolic space [32I . In the case of the vector mode, it is already not obvious 
what square-integrable means and it is not clear which wavelengths we should include in 
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what situation. The results here are of interest regarding this issue because homogeneous 
models seem to correspond to neither k = a/2 nor k = 0, but k = 1. 
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Appendix: Class B metrics in automorphism variables 

For reference, we give a set of metric variables, proposed by Siklos [27[ and their relation 
to the orthonomal frame. First of all, we introduce canonical basis vectors {E^} defined by 

El = dx E2 = (cos ekxdy + e sin ekxdz) 

E3 = e^' (— sin ekxdy + cosekxdz) 

Here the parameter e is taken to be 1 for VII/i, for V and i for VI/j. In this coordinate 
system, the metric tensor can be written in the following form: 



where 



/.2 





e^^' 
e-^A' 

( 

A 



\ 



( 



1 

COS e0 e sin et; 
sine(/) cosed. 



\ 



\ 



1 

e-^ 




B 



1 
1 



$, A and B are called automorhpism transformations since if {^a}a=i,2,?, satisfy the class B 
structure equations , X^/ = (y4$i?)^?Ea also satisfy the same commutation relations. 
By construction, 

ei = aXi, e2 = e'^Xs, 63 = e'^'Xa 
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are orthonormal. Introducing reduced variables 



-e^ ^ ^62 cos e(j) + 636 sin , = 2^^'^'^ (^^^ ^ sin e0 — 63 cos , 



^ 2 

the Ricci tensor in the orthonormal frame is given as follows: 



• 2 

= - - 2^ - 2a\Bl + Bl) - 2A - 2A2 - 2/1^ - 4>\cosh4ii - e'), 
d a 



Roi — 2a + 2aX — e fca0(cosh4/x — e ), 

Rq2 = -2>a^B2 - e^ka^B^e^^", R03 = -Sa^Bs + ka^B2e-^'', 

R,, ^-^ + 2^ + 2a\Bl + Bl) - 2a^ - 2^A - k^a^{coshif, - e^), 
a a'' a 

R12 = {aB2) + 0^2 (3A -ft)- e^aB^^e^^" Ris = (0^3) + 0^3 (3A + /i) + aB2<j)e-^'', 
i?22 = A - /i + (A - /t) (2X - 2a^ - 20^5^ + 2(0^ - k^a^) cosh2/xsinh2/x, 

i?33 = A + /i + (A + |^2A - - 2a^ - 2a^£| - 2(0^ - fc^a^) cosh2/xsinh2/i, 

i?23 = (^-0 + ''''' '^^""'^^ - 2a^B2Bs - e^ka"" - X<i>)(e^^ - e'^e'^n - 2fi^{e^i' + e^e-^'*). 
The LRS speciahsation in section 5 is obtained by setting cf) — /j, — 0. 
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